Detection of the Natural Alpha Decay of Tungsten 
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The natural Q-decay of ^^''W has been unambiguously detected for the first time. The a peak is 
found in a (7, {3 and neutron)-free background spectrum. This has been achieved by the simultaneous 
measurement of phonon and light signals with the CRESST cryogenic detectors. A half-life of 
Ti/2 = (1-8 ± 0.2) X 10^® y and an energy release of Q = (2516.4 ± 1.1 (stat.) ± 1.2 (sys.)) keV 
have been measured. New limits are also set on the half-lives of the other naturally occurring 
tungsten isotopes. 



PACS numbers: 23.60.-|-e, 07.20.Mc, 29.40.Mc 



I. INTRODUCTION 

The a-decay of the naturally occurring isotopes of 
tungsten (W) has been the subject of experimental search 
for many decades. Tungsten is an interesting element 
because a-decay is energetically allowed for all five nat- 
urally occurring isotopes. Mass excess measurements 
show that the available decay energy Q for all these iso- 
topes is low (< 3 MeV) and that the Q values lie in the 
same energy range as /3 and 7 decay due to natural chains 
(see TableCJ . Therefore background suppression and pos- 
sibly event-by-event discrimination against 7's and /3's, 
which may obscure the a signal, are crucial issues for the 
detection of such rare events. 

While I82^,i83w,i84w, and ^^^W are expected to 
have half-lives well above 10^^ y, the isotope ^^°W is ex- 



Isotope 


Abundance [%] 
lUPAC f3] Ref." f4] 


Q [keV] 
Ref. \1] Rcf. \2] 


182W 
183^ 

184W 
186W 


0.12(1) 0.1198(2) 
26.50(16) 26.4985(49) 
14.31(4) 14.3136(6) 
30.64(2) 30.6422(13) 
28.43(19) 28.4259(62) 


2516(5) 2508(4) 
1774(3) 1771.8(2.2) 
1682(3) 1680.0(2.2) 
1659(3) 1656.2(2.2) 
1123(7) 1124(7) 



"Considered by the lUPAC (International Union of Pure and Ap- 
plied Chemistry) to be the best single measurements. 

TABLE I: Representative Isotopic Composition according to 
the lUPAC, best single measurements on isotopic abundance 
from individual references, and decay energies for the isotopes 
of natural tungsten. 
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pected to have a half-life near 10^^ y, and is close to 
observability with present techniques. Most recently, the 
Kiev-Firenze Collaboration, operating CdW04 scintilla- 
tion detectors in the Solotvina mine reported a "first in- 
dication of a possible a-decay of i®°W" with a half-life 
of Ti/2 = Lltllistat.) ±Q.i{sys.) x IQi* years 0. It is 
well known |^ that the scintillation yield for a particles 
is lower than that for /3 or 7 particles of the same energy. 
After pulse shape analysis, an a-peak at ~ 300 keV was 
interpreted as ^*°W a-decay. Due to the low energy res- 
olution (FWHM =110 keV) it was impossible to exclude 
alternative explanations for the peak, hence the result 
was treated as an indication. They also set a 90% C.L. 
lower limit of T1/2 > 0.7 x 10^^ years. Also the ROSE- 
BUD collaboration, operating a 54 g CaW04 scintillating 
bolometer at the Canfranc Underground Laboratory, ob- 
tained a limit of T1/2 > 1.7 x lO" years at 90% C.L 



The high sensitivity and excellent energy resolution of 
low temperature detectors (see for example ^3 ) ^-^^d 
the great efforts made to suppress and understand ra- 
dioactive background in dark matter searches and double 
beta decay experiments make such apparatus sensitive to 
rare nuclear decays. For example, the a-decay of ^^^Bi 
(considered as the heaviest stable isotope) has been re- 
cently detected by a French group developing low temper- 
atures bolometers for Dark Matter direct detection 
CRESST (Cryogenic Rare Event Search with Supercon- 
ducting Thermometers) is a low background cryogenic fa- 
cility primarily devoted to the direct detection of WIMP 
dark matter particles via their scattering by nuclei. Such 
nuclear recoils resulting from WIMP interaction can be 
discriminated from electron background (caused by pho- 
tons or electrons) by measuring phonons and scintillation 
light simultaneously. The energy detected via phonons in 
a cryogenic detector is, to first order, independent of the 
nature of the particle. There is, however, a significant 
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Crystal D: Energy [keV] 

FIG. 1: Pulse amplitude in the light channel versus energy in the CaW04 crystal D from run 28. Three clearly separated 
event populations appear. The solid lines border the three different bands: the upper band is due to 7's and /3's, the lowest 
one to neutrons and the middle one is due to a's. Each peak in the a band has been identified (see text) and is here labelled 
accordingly. 



difference in scintillation yield for nuclei and electrons of 
the same energy. The fraction of energy released in the 
nuclear interaction channel, negligible for photons and 
electrons, becomes important for a's and heavy ions and 
is dominant for WIMPs and neutrons. As a result WIMP 
and neutron induced recoils give considerably less scin- 
tillation light than electrons of the same energy while a 
particles can be clearly discriminated since they interact 
partly with both, nuclei and electrons. This results in 
a pure a spectrum, i.e. one without contributions from 
/3, 7 and neutron events as in the middle band of Fig. 
^ Finally we remark that the energy measured following 
an internal a-decay in a cryogenic detector corresponds 
to the sum of the energies of the a particle and of the 
recoiling nucleus, i.e. the total decay energy Q. 

In this paper we present clear evidence for the a-decay 
of isow from early runs of the CRESST II dark matter 
detectors, originating from the tungsten in the CaW04 
crystals used as the dark matter target. 



II. EXPERIMENTAL SET-UP 

CRESST is an ultra-low background facility operating 
at the Gran Sasso underground laboratories. Detailed 
descriptions can be found in i i. CRESST has de- 
veloped very sensitive cryogenic detectors consisting of 
a dielectric target crystal with a superconducting phase 
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FIG. 2: Schematic of a detector module. It consists of a 
scintillating 300 g CaW04 crystal (phonon channel) and a Si 
wafer (light channel), both read out by a SPT. The set-up is 
surrounded by a reflective foil 



transition thermometer (SPT) evaporated onto one sur- 
face. Particle interactions in a cryogenic detector pro- 
duce phonons which propagate to the surface, where they 
heat the electron system in the SPTs. These thermome- 
ters, biased in the middle of their superconducting-to- 
normal transition near a temperature of 10 milli-kelvin, 
transform the temperature rise of their electrons into a 
relatively large increase of the resistance of the film. The 
resistance change of the low impedance thermometer is 
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then measured via a SQUID. The resistance increase is a 
measure for the energy deposited. 

The set-up for the simultaneous detection of scintil- 
lation light and phonons is described in JP| and shown 
schematically in Fig. |21for one absorber module. It con- 
sists of two independent detectors, each one with its own 
SPT with SQUID readout. The main detector consists of 
a ^ 300g cylindrical CaW04 crystal (0 = 40 mm, h = AO 
mm). CaW04 is a well known scintillator with high light 
yield and very broad emission spectrum peaked at 420 
nm. It is characterized by a high index of refraction (n — 
1.92) and a decay time of 17 /iS at 77 K. The scintillation 
light produced in each target crystal is detected via an 
associated calorimeter consisting of a silicon wafer of (30 
x30 X 0.45)mm'^ volume with a 20 nm thick Si02 layer 
on both surfaces. In CaW04 only a few percent of the 
absorbed energy is transformed into light. To minimise 
light losses the whole module is thus enclosed in a highly 
reflective foil. 

In the data to be presented here, three different 
CaW04 crystals (referred to as crystal B, D and E) were 
used in the Gran Sasso set-up for a total of four different 
CRESST runs (22, 23, 27, 28). For detector operation, 
the temperature of the thermometer is controlled by a 
dedicated heater. Additionally the heater is used to in- 
ject test pulses which monitor the long term stability of 
the detector 0. In particular, in the runs 27 and 28 a 
large voltage pulse was produced, reaching the saturation 
region of the superconducting transition specifically to 
monitor the detector response above 1 MeV. This pulse 
was sent, at a rate of 0.5 Hz, throughout the measur- 
ing period. A record length of 4096 samples and a time 
base of 40 /zs, resulting in a time window of 164 ms for 
recording an event, was chosen for all the runs presented. 



III. OFF LINE ANALYSIS 

CRESST detectors are optimised to be most sensitive 
in the energy region relevant for dark matter direct de- 
tection (<200 keV). In this energy range the detector 
response is mapped out with electric heater test pulses, 
and the 122 keV line from a ^^Co source provides the 
calibration of the test pulses in terms of a 7 equivalent 
energy @. A pulse shape template is obtained by aver- 
aging pulses from the 122 keV line. This template is used 
in a fit to the actual signal pulses in order to accurately 
determine the signal amplitude. This is valid as long as 
the detector response is in a linear regime. If the re- 
sponse becomes non-linear the pulse shape changes and 
the same template cannot be used over the full energy 
range. 

Typically, the pulse shape for the CaW04 detectors 
can be described by Trise = 1-1 ms and Tdecay — 30 ms. 
However, the temperature rise caused by a-decays is usu- 
ally beyond the dynamic range of the SPT and for precise 
energy determination the signal duration rather than the 
pulse height must be used. Therefore two new analysis 



methods, described in [13, were developed in order to 
reconstruct the higher energies (MeV region), for which 
the detector response is non-linear. 

In the first method (Method A) a template from pulses 
in the linear region (for example the 122 keV template 
from the calibration) is used to fit all measured pulses 
with an amplitude below the voltage at which the re- 
sponse becomes non-linear. For all pulses with an ampli- 
tude exceeding that voltage the fit is truncated and the 
upper part is reconstructed from the template. This is 
shown in Fig. 13 The resulting amplitude spectrum is 
then linear by construction. 

If the pulse duration is too long the information from 
the decaying part of the pulse is outside the recording 
time window and the method fails. To overcome this lim- 
itation a second method was developed (Method B). The 
template from the linear region is fitted with a function 
following the model described in This assumes an 

exponential rise and two exponential decay components. 
Once the parameters of the model function are deter- 
mined for the linear region, a mathematical description 
of the changes in pulse shape in the non-linear range is 
needed. At higher energies the pulse shape varies accord- 
ing to the non-linearities of the superconducting phase 
transition. These changes can be described via successive 
approximation to the model function. The higher order 
correction coefficients are extracted by fit ting templates 
of different energies up to 2.31 MeV (see [13), enabling 
the energy spectrum to be reconstructed up to the high- 
est recorded a-signal [a — a cascade at ~15 MeV). 



IV. BACKGROUND 

Due to the powerful discrimination technique CRESST 
II detectors offer, contaminations from natural decay 
chains and other a-unstable isotopes have been identified 
by their a-decays with a sensitivity of ^ 1/iBq/kg. The 
results of this analysis have shown very good agreement 
with the results obtained using different techniques com- 
monly employed to determine the presence of impurities 
in crystals. These include ICPMS"'^ (Inductively Coupled 
Plasma Mass Spectrometry), HPGe^ (High Purity Ger- 
manium 7-spectroscopy) and X-ray luminescence'^ tech- 
niques. 

Evidence for contaminations due to a-unstable rare 
earth elements has been found in all the measured detec- 
tors. Rare earth doped crystals are commonly produced 
for laser applications, therefore these impurities are likely 
to have been introduced during the production of the 
crystals. In particular an a peak at 2.31 MeV has been 
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FIG. 3: Fit of two events from crystal D (run 28) with a template obtained by averaging pulses from the 122 keV peak. The 
fit is truncated at 0.7 V. For pulses of amplitude below this value the fit procedure is standard (left plot). For pulses with 
an amplitude exceeding this value, the higher part of the pulse shape (above the horizontal dotted line in the right plot) is 
distorted due to saturation. Therefore the fit is truncated, i.e. fit only to data points below 0.7 V, and the amplitude of the 
pulse is reconstructed from the template. 





Grystal B 


Crystal E 




semi- 






semi- 




quantitative 


quantitative 


quantitative 




analysis 


analysis 


analysis 


Natural 


Sample I 


Sample II 


Sample I 


Sample II 




Element 


[ppb] 


[ppb] 


[ppb] 


[ppb] 


[ppb] 


Nd 


555 ± 167 


773 ± 232 


683 ± 137 


994 ± 199 


1100 ± 330 


Sm 


<22 


11 ± 3.3 


9 ± 2 


13 ± 3 


<6 


Eu 


<8 


4 ± 1.2 


4 ± 1 


5 ± 1 


<2 


Gd 


2898 ± 869 


4084 ± 1225 






1000 ± 300 


Er 


1328 ± 398 


1877 ± 566 






940 ± 283 


Dy 


<15 


<3 






<3 


Hf 


<22 


< 4 






<3 


Os 


<29 


< 5 






<4 



TABLE II: Results from the ICPMS analysis at the Gran Sasso facility on rare earths, hafnium, and osmium. For the 
semiquantitative analysis the calibration curve is determined with a multi-element standard which may be different from the 
one to be analysed. On some elements, in particular on samarium, a quantitative analysis was performed. In this cases the 
calibration is obtained using known standards of the element to be analysed. All the samples are dissolved in a HNO3 solution 
in the microwave oven with different exposures for samples I and II. 



detected and identified as the a-decay of ^^^Sm. This 
isotope is naturally occurring with 15% isotopic abun- 
dance and is known to be a pure a unstable radionuclide 
with Ti/2 = 1-06 X 10^^ years. The intensity of this peak 
varies for the different detectors, implying different lev- 
els of contamination in each crystal. Crystal B showed a 
count rate of (1.7 ± 0.1) counts/h which translates into 
(13 ± 1) ppb of natural samarium in the crystal. The 
rate on the same peak measured with crystal E corre- 
sponds to a contamination of (5.5 ± 0.4) ppb of natural 
samarium. 

ICPMS analyses were performed in the Gran Sasso fa- 
cility on crystals B and E giving results, shown in Table 
nil in agreement with our analysis of the a-peak intensi- 
ties. Evidence of ^''^Nd (Q = 1905.1 keV, isotopic abun- 
dance (i.a.) = 23.8%, Ti/2 = 2.29 x lO^^ y) and ^^'^Gd 
(Q = 2205 keV, i.a. = 0.2%, T1/2 = 1.08 x lO^'* y) [H 



a-decay has also been observed and the results are in ac- 
cordance with the ICPMS analysis. Furthermore, X-ray 
luminescence analysis on crystal B showed the presence 
of the characteristic emission peak of Gd and Er |3| , the 
two most abundant rare earth contaminants (see Table 
III. 



Finally, alternative candidates such as ^^*Hf (Q =2496 
keV, i.a. = 0.162%, T1/2 = 2.0 x lO^^ y) and ^^^Os 
(Q = 2822 keV, i.a. = 1.58%, T1/2 = 2.0 x lO^^ y)[l3 
were considered to cause the 2.31 MeV peak. They were 
discarded, however, due to the unreasonably high atomic 
concentration needed to produce the detected rate. 
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FIG. 4: Top: Fitted peak position versus known energy. The energy scale for the y-axis is fixed by the Sm peak while the 
position of the ^^^Th and ^'^^U peaks results from a simple linear extrapolation of the calibration for both methods A and B. 
Bottom: Residuals of the fitted position from the expected position. For method A, the ^•^^Th peak results at (4081.3 ± 1.7) 
keV (AE = (28 ± 4) keV (FWHM)) and the ^^^U at (4269.4 ± 2.2) keV. For method B, the peaks lie at (4076.3 ± 0.8) keV 
(AE = (14.1 ± 1.7) keV (FWHM)) and (4265.9 ± 1.1) keV respectively. The calibration of the amplitude spectrum resuhing 
from method A is much more accurate; however, the resolution is less compared with method B. 



V. CALIBRATION 

Once unambiguously identified, ^^^Sm was used to cal- 
ibrate the full a spectrum. This allowed a precise iden- 
tification of the U-Th peaks, which were then used to 
determine the accuracy of the calibration for both re- 
construction methods A and B. In Fig. 2|thc linearity 
of the response functions for method A (top left figure) 
and method B (top right figure) is shown for the data of 
run 28. On the x-axis the central values of the tabulated 
energies Q of the ^-^^Sm (Q= 2310.5 ± 1.1 keV), "^^^Th 
(Q= 4082.8 ± 1.4 keV) and "^^^V (Q= 4270 ± 3 keV) 
peaks are given. On the y-axis their fitted peak posi- 
tion is plotted. The energy scale for the y-axis is fixed 
on the ^^^Sm peak at 2310.5 keV while the position of 
the ^^^Th and ^ss^j peaks results from a simple linear 
extrapolation of the calibration. At the bottom of the 
figure the residuals of the fitted position with respect to 
the expected position are shown. The calibration of the 
amplitude spectrum resulting from method A is much 
more accurate, but less information can be used to re- 
construct the signal amplitude compared with method 
B. This results in a slightly worse energy resolution AE 
= (12.9 ± 0.3) keV (FWHM) on the "^Sm peak, and 
it fails to completely reconstruct the higher part of the 
U-Th a-spectrum. Method B, on the other hand, us- 
ing a larger amount of information from the pulses than 



method A, gives a much better resolution (AE = (6.7 
± 0.1) keV (FWHM) on the ^-^^Sm peak), but the ap- 
proximation of the pulse shape description introduces a 
higher systematic error on the linear extrapolation of the 
calibration. In the following, method B has been used 
to fully reconstruct and identify the entire a spectrum, 
while method A was applied as an independent check on 
the peak position of the ^*°W evidence. In particular, 
method A is used to extract the measured energy for the 
i*°W Q-value. 



VI. RESULTS 

Crystal D showed the best performance as a detector 
(energy resolution: AE = 1.8 keV (FWHM) for 122 keV 
gammas, and AE =(6.7 ± 0.1) keV for 2.31 MeV alphas, 
using method B) and it has the longest exposure (11.745 
kg days in run 27 and 12.268 kg days in run 28). In run 28 
the performance of the module was improved by replacing 
the Ag reflective foil used in run 27, which showed a 
^^°Pb contamination, by a highly reflective (99% at 420 
nm) scintillating polymeric foil. It is these data which 
we discuss in full detail. A peak was observed at 2512.8 
± 0.4 (stat.) ± 4.3 (sys.) keV with a resolution AE 
= (6.5 ± 0.7) keV (FWHM). This is shown in Fig. 
The position in energy obtained by reconstructing the 
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FIG. 5: Alpha events measured in Run 28. In the top figure the energy spectrum between 1.7 and 2.7 MeV is shown with 3 
keV bins. On the vertical axis the number of counts are displayed on a logarithmic scale. In the bottom figure the fit to the 
^*°W Qf-decay is shown. No background counts are present in a wide energy interval around the peak. 



spectrum with method A is (2516.5 ± 1.4) keV. A half-life 
Ti/2 = (1.7±0.2) X 10^^ y is obtained. This is consistent 
with the value from run 27 on the same crystal and on 
crystals B and E, as shown in Table Hill The consistency 
of the results from different crystals tends to confirm that 
the signal originates from the ^^^W decay and not from 
unknown impurities which may be expected to vary for 
each individual crystal, as seen in Section 4. Furthermore 
ICPMS measurements performed on these two crystals 
reported in Table ITTl exclude alternative explanations for 
the tungsten peak, such as -'^^'^Hf. 

Finally, the four measurements were added, resulting 
in a spectrum for a total exposure of 28.62 kg days. The 
data from each run have been analysed with method A, 
calibrated individually and then summed. This results 
in half-life of T1/2 = (1.8 ± 0.2) x IQ^^y, which is con- 
sistent with the previously published limits. Including 
the uncertainty on the position of the ^^^Sm line, a peak 
of energy (2516.4 ± 1.1 (stat.) ± 1.2 (sys.)) keV is 
found with a resolution AE = (18 ± 2) keV (FWHM), 
as shown in Fig. |S1 This result is consistent with expec- 
tations based on the mass difference of ^^"W and ^^^Hf 
from Ref.T'l, but lies 2cr away from the latest update ^ 
(see Table (J) . We note however that most of the input 
data used to evaluate the atomic masses consist of rela- 
tive measurements which set a relation in mass or energy 
among two or more nuclei. In the diagram showing the 
relations among input data (Fig 1 of Ref.Q]), ^*°W and 
"'^^^Hf are only weakly connected and an underestimation 
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FIG. 6: Sum of the four different runs, each analysed inde- 
pendently with method A. The resulting peak is fitted giving 
an energy of (2516.4 ± 1.1) keV with resolution FWHM = 
(18 ± 2) keV. Few background events due to a surface con- 
tamination of the Ag foil used in run 27 are visible. The 
scintillating foil used in runs 22, 23 and 28 efficiently vetoes 
all these events (see text). 



of the errors is plausible. Our result represents the miss- 
ing direct connection between ^^"W and ^^^Hf. 

Limits on the half-lives for the non-observed a-decays 
of the other four naturally occurring tungsten isotopes 
were also calculated. For this purpose the data from run 
27 could not be used. In that run, a Ag foil was used to 
surround the crystal, while in the other three runs the 
scintillating polymeric foil was used. Any contamination 
of a-unstable isotopes, such as ^^"Po on the surface of the 
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Crystal Run 


Exposure Number 


Half-life 


number 


Wff Havsl 


of counts 


M 




D 28 


12.268 


35± 4.3 


Tl/2 


= (1.7 ±0.2) X 10'** 


D 27 


11.745 


28.5 ± 5.6 


Tl/2 


= (1.9 ±0.4) X 10^® 


E 23 


3.467 


9 


Tl/2 


> 1.1 X 10'** (90% C.L.) 


B 22 


1.14 


4 


Tl/2 


> 6.8 X 10'^ (90% C.L.) 


All Crystals 


28.62 


75.5±8.7 


Tl/2 


= (1.8 ± 0.2) X 10'** 



TABLE IIL Half-life for the a-decay of W obtained from the present work. The results are given independently for the four 
different runs and for the total exposure. For the runs 27, 28, and the total exposure, the number of counts is determined by 
a Gaussian fit to the peak. Errors are given at la. For run 22 and 23, due to the small number of counts, the total number of 
events contained within ±3(t of the expected energy position is given. 



Isotope Half-life [y] 

This work Previous Previous ]Jj 

'"'^W Tl/2 > 7.7 X 10^' > 1.7 X 10^" > 2.5 x 10'' 
i83w Tl/2 > 4.1 X 10^' > 0.8 X 10^° > 1.3 x lO'" 
i«4W Tl/2 > 8.9 X 10^' > 1.8 X 10^° > 2.9 x lO" 
i86w Tl/2 > 8.2 X 10^' > 1.7 X 10^" > 2.7 x lO'" 



TABLE IV: Lower limits (90% C.L.) on the half-lives of the 
non-observed a-decays of the other four naturally occurring 
tungsten isotopes, obtained with a total exposure of 16.875 
kg days. For comparison the best previously published limits 
are also shown . 



foil or on the crystal may results in a events of degraded 
energy from the U-Th chain energy, going down to ~100 
keV. These were observed in run 27 with a rate of 0.087 
counts/lOOkeV/kg/day both between ~500 keV and the 
Nd peak at 1.9 MeV, and between the '8°W and the 
^^^Th peaks. The polymeric foil used in all the other runs 
did not show any surface contamination. Furthermore if 
the scintillating side of the foil faces directly the crystal, 
as in run 28, it efficiently vetoes all the surface events 
due to the light emitted by the polymeric foil itself. The 
limits obtained from the sum of runs 22, 23 and 28 are 
reported in Table Hvl 



VII. CONCLUSIONS 

We have observed the natural a-decay of '*°W with 
a half-life of T1/2 = (1.8 ± 0.2) x lO^^y and a Q-value 
of (2516.4 ± 1.1 (stat.) ± 1.2 (sys.)) keV. In addition, 
the lack of any signal from the other tungsten isotopes 
sets new limits on their half-lives. These results offer an 
improvement of roughly a factor of 50 over the best previ- 
ously published limits. Finally we would like to note the 
power and potential of this cryogenic technique, which 
here enabled us to achieve a very high sensitivity in an 
application to nuclear physics. 
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